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Abstract 

The electromagnetic field can be expressed in terms of two complex potentials a, , which are 
related to the Debye potentials. The evolution equations for a, f3 are derived, which are separable 
either in parabolic coordinates (leading to the radiation fields) or in radial coordinates (multipole 
fields). Potentials corresponding to focused wave fields as well as plane waves are discussed. A 
conserved radiation density can be constructed in terms of these potentials, which is positive 
(negative) for positive (negative) helicity radiation. 

PACS nos. 4110 classical electromagnetism 

4110H electromagnetic waves - theory 



I. INTRODUCTION 



The source- free electromagnetic field has only two degrees of freedom per space-time point. Some 
economy is achieved by reducing the six components of the electric and magnetic fields (E, B) to 
the usual 4 potentials ($, A) which satisfy the covariant Lorentz condition d t Q> + V ■ A = , or the 
<E> component may be dispensed with provided that the three potentials A satisfy the non-covariant 
Coulomb gauge V ■ A = . The most economical way to express the two degrees of freedom of the 
free electromagnetic field is in terms of the two real Debye potentials 1,2 (tjji, ^2) ■ Here we introduce 
a pair of complex potentials {a, (3} from which the electromasgnetic field can be derived via (8) 
below. Some advantages of using these potentials are: a conserved radiation density representing 
the difference of positive and negative helicity radiation can be constructed, they accomodate 
singularities in the field (charged particles), and interesting focused wave solutions arise naturally 
when solving the evolution equations (equation (A7)below). 

We first review the Debye potentials to which {a,/3} are related. The electric and magnetic 
fields are expressed in terms of the Debye potentials by 
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where ip are solutions of the wave equation 



V 2 ) ^ = 0. 



(2) 



The standard text-books may not refer to the Debye potentials by that name 3 , but use either ip 
or L 2 ip in their analysis of multipole radiation fields, where tp are spherical solutions of the wave 
equation - see for example pp. 432-433 of Reitz 4 , or pp. 745- 746 of Jackson 5 . The paper by Boukamp 
and Casimir 6 discusses of the various essentially equivalent approaches to multipole radiation, while 
Refs. 2,7 show that any field (E, B) outside the source region can be expressed in terms of the 
potentials %px , ip 2 ■ 



We express (1) more concisely by defining the complex field F = E + iB , then 
F = [-(x x V)d t + iV x (x x V)] + #2) 

F= [idt + Vx] L4> (3) 

where L = — ir x V is the angular momentum operator and ip = —(ipi + ifa) ■ Conversely 

r-F = r- VxL^ = iL 2 ip (4) 

and the operator L 2 is, in principle, invertible. 8 The source-free Maxwell's equations are 

[id t -Vx]F = V-F = 0. (5a, 6) 

It is readily verified that substituting (3) into the above yields an identity, given (2). 

Despite the fact that the two Debye potentials (or the single complex Debye potential ip ) 
precisely contain the two degrees of freedom of the electromagnetic field, they are seldom used 
outside magnetostatics and multipole field analysis. Some disadvantages spring to mind: 

(a) the field F = E+iB is expressed by second order differential operators acting on the potentials, 
so then for example the usual energy density of the field (F* F) is a cumbersome expression in 
terms of tp , 

(b) there is no regular closed form potential ip representing a plane wave field. The standard 
texts 9 express the plane wave potential as an infinite sum over I of the spherical harmonics Yi,±i , 
although this actually represents plane waves outgoing in each direction from the (9 = tt/2) plane, 

(c) the Debye potential formalism is ill suited to coping with point singularities representing charged 
particles in the field, i.e. there is no regular potential %p representing the monopole field at the 
origin, as can be seen immediately from (4). See Ref.f for irregular potentials ip representing the 
monopole field. 

We will now show that, given (3), the field F can be expressed in terms of two complex potentials 
a, (3 acted on by first order differential operators. These potentials (a, (3) admit charged particle 
solutions, for static and moving charges, and the potentials representing the plane wave fields are 
fairly simple. 

From (3) we calculate the field (F + ixx F) in terms of ip , where x = * . We use the operator 
identities 2 

V x L = (-ixV 2 + iV(d r r) ) , x x L = i (Vr - <9 r x) , x x (V x L) = -L(<9 r r) . 

Then 

F + zx x F = [i d t L - ixV 2 + iV(d r r)} i/> + i [ -(Vr - d r x)d t - L(0 r r)] V 

F + zx x F = [V + xd t - -L] a (6a) 

r 

where 

a — i(d r rip — d t rip) . (7a) 

To derive (6a) we have used (2). Similarly 

F - ix x F = [V - xd t + -L] 13 (6b) 

r 

with 
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(3 = i{d r rip + d t rip) . {7b) 

Note that given (F + ix x F) there is no straightforward algebraic manipulation to derive F due 
to the fact that the matrix (/ + ixx) is singular. We add (6a,b) obtaining 

2F = [V + m - -L] a + [V - xd t + -L] (3 , (8) 

r r 

which is the basic formula expressing F in terms of {a, (3} . Dotting (6a,b) with x yields the 
constraint 

(3 r + d t )a = F r = (3 r - dt)0 (9) 

where F r = x • F . The operators (d r + dt) , (d r — dt) are the radial differential operators on the 
future, past lightcones respectively. Given the radial field F r (t, x) , then (9) enables us to calculate 
a, (3 as follows. To carry out the integration on the future (past) lightcone we introduce the 
advanced (retarded) times S and T defined by 

S = t-r, T = t + r. (10) 

The hypersurfaces S = So, Si... are future lightcones centered at the origin, while the hypersurfaces 
T = Tq,Ti... are past lightcones. Then F r (S) or F r (T) are obtained by substituting t = S + r or 
t — T — r into the expression for F r (t) . A solution of (9) is then 

a(S)= [ drF r (S) P(T)= I drF r {T). (11) 
Jo Jo 

Maxwell's equations can alternatively be formulated from the start using advanced or retarded 
time, 10 then a or (3 is found from F r by a straightforward radial integration. 

We now find the evolution equations for a, (3 , noting that a = i(d r r—dtr)ip does not satisfy the 
wave equation like ip , because the operator (d r r — d t r) does not commute with the d'Alembertian 
operator □ = (<9 t 2 — V 2 ) . The d'Alembertian may be expressed as follows, 



(d 2 -v 2 )v = i-kd r + d t )(d r -d t )r + 1l 2 ]^ = o. (12) 
3 



Now multiply (11) by ^(d r — d t )r obtaining 

[~(d r - d t )r 2 (d r + dt) + ^L 2 ]i{d r - d t )r1> 

= [9 2 -V 2 - \]a = 0. (13a) 

so that the evolution operator for a is a modified d'Alembertian, having the extra term — ^dt ■ 
Similarly 

[3 2 -V 2 + -$]/? = 0. (136) 

We see from (8) that adding any function f(t — r) to a , or adding any function g(t + r) to 
(3 leaves F unchanged. Also the evolution operators (13) for a, (3, annihilate f(t — r),g(t + r) 
respectively. So the gauge transformations for a, (3 are 

a — > a' , a' = a + fit — r) 

(14) 

(3' =a + g(t + r). 
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II. THE POTENTIALS {a, (3} CORRESPONDING TO RADIATION 
FIELDS 

We first establish some symmetries between a, (3 . Consider the parity transformation V : x — > — x 
and the time reversal transformation T : t — * — t . Under the combined transformatioon VT the 
field transforms as 11 

PTF(t,x) = _F(-t,-x). 

For brevity we will write a bar over any variable signifying its VT transform, e.g. F = VT F , a = 
VT a etc. Both a, (3 are scalars under VT , so that a(t, x) = a(— t, — x) . Maxwell's equations are 
invariant under the VT transformation, so that if F satisfies Maxwell's equations then so does 
F' = F . In some cases of interest the field is even or odd under VT , by which we mean F ev = F ev 
and F dd — —F dd (from now on we will use the terms even and odd meaning even or odd under 
VT .) Recall the fields (F + ixxF), (F — zx x F) which are derived from a, (3 respectively. Then 
if the field is even/odd, 

(F ev -ixxF«) = ±VT (F +ixxF ev ) . (15) 

v odd odd ' x odd odd ' v ' 

In either of these cases the field is effectively given by one potential, a say, because after calculating 
(F + ixx F) from a , one can write down (F — ixxF) using (15). Applying the VT transformation 
to equations (6): 

F - zx x F = -[V - Std t + -LI a F + ix x F = -[V + x<9 t - -Ll (3 , 

r r 

so that 

F — > F is equivalent to {a,f3} —> {—(3, —a} . (16) 

From (16) we see that 

{F ev : (3 = -a} {F odd : /? = +«}. (17) 

We now find the potentials {a, f3} corresponding to the radiation fields. Equations (11) 
enable us to calculate the potentials for a known field - which approach we follow in this section. 
Alternatively we can look for solutions of the evolution equations for the potentials (13), satisfying 
the constraint (9) - which approach we take in appendix A. It turns out that these evolution 
equations (13) are separable in either spherical or parabolic coordinates - the spherical solutions 
yield the potentials corresponding to the multipole fields (which are singular at the origin) , while 
the parabolic solutions are potentials for non-singular fields, including the plane wave fields. 

Consider the plane wave field F propagating in the x% direction, the left circularly polarized 
(positive helicity) field F L is 

F L = {l,i,0} e ikx3 - ikt F rL = (xi + ix 2 ) e ikx? '~ ikt (18) 
and its VT transform the right circularly polarized (negative helicity) field F^ is 
F R = F L = -{1, i, 0} e -(^x 3 -ikt) FrR = + j£ 2 ) e -(ikx 3 -ikt) 
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Then recalling (10) we integrate F r L(S) = (xi + ix?) e lks e lkr O- x z) to find a, and follow the 
equivalent procedure for (3 , obtaining 



X\ + 1X2 \ r ikx 3 _ ikr] £_ 



-ikt 



1 1 j 



. : r:i - ;/ '' ifet (19) 

The second term in the square brackets does not contribute to the field but its absence would make 
the potentials singular due to the factors 1/(1 f x 3 )) in a, /?. Direct calculation verifies that 

[V + x<9t - -L] a L = { 1 + £ 3 , i + ix 3 , -x x - ix 2 } e ikx? -- ikt 

which is Fl +ix x Fl , noting that the operator [V + xdt — ^L] commutes with the ( ^pj^ 2 ) factor 
in ai, , and annihilates the term e lfe ( r -*) . Similarly it can be verified that [V — x.d t + ^L] /?l = 
F L - iSt x F L . 



As Fr = Fl , then recalling (16) 

( 20 ) 

f>H = -a L = (^±^)[e-^-e^] e —. 
The potentials for any plane wave propagating in the £3 direction are a linear superposition 
of {ccljPl} , {ctR,f3fi} . Both {o.l,[3l,} , {aR,f3ji} are eigenstates of L 3 , the angular momentum 
operator in the direction of propagation, with eigenvalue +1 . 

We obtain the potential ai(k) for a left circularly polarized wave propagating in the k = k/|k| 
direction with wavelength 27r/|k| , by first defining the null complex vector h satisfying 

hxk = ih h*h = 2, (21) 
which so defined is unique up to a phase factor. Then 

{a L (k), /9 L (k)} = { - (-^-) [e lk ' x - e l ' k 'l , (-^-) [e lk ' x - e"^ } ^ , (22) 
r — k • x r + k- x > ik 

and {«i?(k), /3i?(k)} = {— /3z,(k), — a.£,(k)} . We check that az,(k) , /?z,(k) are eigenstates of (k • L) 
the angular momentum operator in the direction of propagation, with eigenvalue +1 : 

k L (h x) = k (-ix x h) = -zx • (h x k) = -ix • (ih) = h • x . (23) 

In the next section we investigate the orthonormality of the basis potentials! a (k), /3(k)} and 
Wk'),/3(k')}. 

We next briefly consider the multipole radiation fields and their corresponding potentials, 
noting some symmetry relations under the VT transformation. The Debye potential tp in this 
case is a spherical solution of the wave equation, then the potentials {a, 13} follow from (7). 
Consider 

ipim = 3i(kr) Y lm {6, (j>) cos(fci) (24) 

where ji is the spherical Bessel function of order I , and Y; m are the spherical harmonics. This ipi m 
when substituted into (2) yields F Jm , the magnetic multipole field of order (/, m) . (Multiplying the 
potential by i yields the corresponding electric multpole field iFi m ). Then from (7) the potentials 
{ai m ,(3i m } for the magnetic multipole field of order (l,m) are 

aim = i(d r - d t )[rji(kr) Y lm (6,4>) cos(Art)] 

A™ = i{dr + d t )[rji(kr) Yi m {9, 0) cos(fci)] . 



We see that 

Plm = (-1)' ai m 

following the parity of Yi m . This means, recalling (17), that F; m is odd or even, depending on 
whether I is an even or odd number. (If we had considered the potential ip of (21) with a (sin(fci)) 
instead of a (cos(kt)) factor, then (3' lm — (— 1)' +1 a' lm .) Then, as discussed previously, F im can be 
calculated from either one of the potentials a, (3 . 

III. THE CONSERVED DENSITY p AND THE ENERGY E 
A. The radiation density p . 

The electromagnetic field is the classical "first-quantized" version of the massless spin-1 photon 
field, and as such should have a conserved density corresponding to the well-known densities for 
the Klein-Gordon (spin-zero) or Dirac (spin-|) fields. The only text on classical (non-quantized) 
electromagnetism which discusses a radiation density, or an inner product, for the electromagnetic 
field that I have come across is the book by Good and Nelson, 12 otherwise one has to turn to 
texts on QED such as that by Schweber. 13 For an excellent recent review on this subject, see 
Ref. 14. This inner product for the electromagnetic field involves the non-local (integral) operator 

(v 2 )- 1 / 2 . 

The evolution equations (7) for a , (3 can be expressed in the form 

(V-xQt) • (V+x<9 t )a = 

(25) 

(V + xft) • (V-x«9 t )/3 = 0. 

It follows that 

(V + x5 t ) • [a (V - W t )f3* - a* (V - Zd t )(3] 

- (V - W t ) ■ [13* (V + Zd t )a - (3 (V + x<9 t K] = . (26) 

Collecting terms and multiplying by (— i) we can write (4.3) in the form of a continuity equation 
d t p + V • J = , where 

P = ~[a (d r - d t )(3* - a* (d r - 8 t )f3 + f3*(d r + d t )a -f3(d r + d t )a*] 

p = -'-{a - (3) F; + CC =- l -[{a-(3)(d r + d t )a*} + CC (27) 

J = - l - [a (V - W t )f3* - a* (V - Hd t )(3 - (3* (V + ±d t )a + (3 ( V + W t )a*} (28) 

where we have used (9), and CC stands for the complex conjugate terms. Let us substitute the 
{a L ,(3 L } of (19) into (27): 

p L = - i- ( - ( ^ 1 ± ^ 2 ) [e^ - e* kr ] - + li 2 ) [e^ - e~ ikr ] ) ((x, - ix 2 ) e^) + CC 

= ((1 + x 3 ) [1 - cos(fcr - kx 3 )j + (1 - x 3 ) [1 - cos(£x + kx 3 )} ) > . 

Inserting {o.r,(3r} = — {/3_r,Q!_r} into (27), we obtain pr = —pL . Hence pl is non-negative and 
Pr is non-positive. 
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The density p has close similarity with the Klein-Gordon density for a spin-zero particle: 
Pkg = ~^{ ( t>*^t4> ~ (t>dt4>*) ■ I n this case solutions <f> with time dependence e~ zujt with to positive 
are regarded as particle solutions, those with time dependence e +lujt are regarded as anti-particle 
solutions. Analagously we can regard the right circularly polarized (negative helicity) field as the 
'anti-photons' of the left circularly polarized (positive helicity) photon field. There appears to 
be no comparable radiation density in the literature, although of course one can project out the 
different polarizations when the field is expressed as a Fourier integral. 15 ' 16 



We now investigate the orthogonality of the potentials for two different fields Fi and _F 2 . It 

Pi 



will be convenient to represent the potentials for the field Fi by the 2x1 matrix U\ 
then we construct the following indefinite scalar product space 

(*7i|t/ 2 ) = J d 3 x Pl2 

where pi2 is the conserved density 

if? (idt + <r 3 d r )u; - uT (id t + a 3 d r )[/i 



(29) 



Pl2 = -4 



Pl2 = ~l [(«1 " Pi) K2 ~ K " 0*2) Frl] ■ 







(30) 

Note that when 



Here U T = [a 0\ is the transpose of U , and <r 3 is the Pauli matrix 
U1 — U2 — U, then P12 = p . In appendix B we show that 

(C/ L (k)|C/ i (k / )) = (27r) 3 fc ( 5(k-k') 

(U R (k)\U R (k')} = -(2tt) 3 k 8(k - k') 

(U L (k)\U R (k')} = . 

The relations (31) enable one to project out the positive/negative helicity states of momentum k, 
and show that p integrated over all space is the amount of left minus right circularly polarized 
radiation. 



(31) 



B. The energy density E 

As a consequence of Maxwell's equations (5) 

S t (VF*) + V-(^FxF*)=0 

where E = \~F ■ F* = i(E 2 + B 2 ) is the usual energy density, and |FxF* = ExB is the 
momentum density. To simplify F • F* in terms of a, (3, we note the following: that from (8) we 
can express F as 17 



F = V T a + V x /3 + xF r 
with F r expressed in terms of a, ft by (9), and the operators Vf, Vj, are defined as 



(32) 



V r = ^(V-xa r +jxxV), Vj. = (V T )* = i (V-x<9 r -jxx V) . (33) 



These operators have the property that for any scalars x, ( 

v T x • v T c = v a • v x c = x • V TX = x • v iX = 

V T x • V X C = i[(x x Vx) • (x x VC) + ix • (Vx x VC)] ■ 



(34) 
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Then substituting in (32) into (F • F*) , and with the aid of (34), we find 

2F = F-F* = [V T a + V x /3 + xF r ] • [V t a* + V T /3* + xF r *] 

2F = |V T a| 2 + (Vi/31 2 + |F r | 2 (35) 

with F r given by (9). We will further expand (35) as we shall see that E contains a divergence 
term, which can be removed such that the resulting energy E' is still conserved. Again using (34) , 

4E = |x x Va\ 2 + |x x V/3| 2 + 2 |F r | 2 + ix • [(Va x Va*) - (V/3 x V/3*)] , (36) 

and the last terms of (36) can be expressed as iV • [—a (x x Va*) + /3 (x x V/3*)] , so that the 
modified energy density E' 

AE' = |x x Va| 2 + |x x V/3| 2 + 2|F r | 2 (37) 
is also a conserved non-negative density, so is also a plausible candidate for the energy density. 



IV. THE POTENTIALS CORRESPONDING TO THE FIELD OF A CHARGED PARTICLE 

First consider the static case, then the constraint (d r + d t )a = (d r — d t )(3 implies j3 — a , then 
from (8) we obtain 

F stc = Va stc 

so that 

Re[a] = Im[a] — $ m 

where $ is the usual scalar potential, and $ m is the magnetic scalar potential. The potentials 
$ , $ m for the electrostatic, magnetostatic fields are well known. The energy E' for the static 
field is E'= i|Va| 2 . 

Magnetic monopole fields are also accomodated: the potential for a magnetic mono-pole at 
position a is just a = (3 = i q m /\ x — a | • Note that if we use (11) to derive a from the radial field 

F r — q(r - x • a)/|x - a| 3 

of a charge at position a, we obtain a = J Q r dr F r — — (g/|x — a|) + (g/|a|) . For non-radiation fields, 
we can lose the constant of integration by instead obtaining a from F r , so that (9) is satisfied, as 
follows 

y(S) = f drF r (S). (38) 

•J oo 

We will use (38) to determine the potentials {a, (3} for a uniformly moving charge. The field of a 
charge passing through the origin in the x\ direction is 18 

F — {xi - vt, x 2 — ivx 3 ,X3 + ivx 2 } F r = (39) 

where r' = [y 2 {x 1 — vt) 2 + x\ + x 2 } 1 / 2 . Substituting t = S + r into the expression for F r , 

rp r q\ _ q~f{r - y Xl S - v Xl r) 

r{ ) [^( Xl r-vS- vr) 2 + (x 2 +x 2 3 )r 2 ] 3 / 2 [ ' 



a 
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and integrating over r we obtain 

<n[r(l - X\v) + S] 



a 



S[y 2 {x\r — vS — vr) 2 + {x\ + a; 2 )?" 2 ] 1 / 



where t' — j(t — vx\) . This potential is non-singular except at the charge position, because when 
t — r, r' — 7(r — vx\) , and ( * t ~^ ) = 7 . The field F of (39) is an even field (which is only the case 
because the particle is passing through the origin), and so from (17) 

(3=-a = -q{ — )—■ 
v t + r r 

V. OUTLOOK 

We have calculated the potentials {a, (3} corresponding to various fields: the simple relation (9) 
between the potentials and the radial field F r makes the calculation of the potential from a given 
field quite straightforward, or one can solve the equations for the potentials from which one then 
calculates the field, as in appendix A. Although {a, (3} are scalars under rotations, they have 
complicated transformation properties under the Lorentz transformations - we will discuss these 
transformation properties elsewhere. (For the Lorentz transformation of the Debye potential ip , 
also complicated, see the paper by Monroe. 19 ) The stationary solutions for the potentials with time 
dependence e~ lkt (e +lfc * ) correspond to the left (right) circularly polarized waves: the orthogonal- 
ity properties (31) suggest an alternative approach to second quantizing the electromagnetic field, 
which we hope to address in the future. 

Under the duality transformation (E — > B , B — > — E) , or F — > iF , the transformation of 
the potentials is just {a, (3} —> i {a, (3} . The absence of any duality transformation for the usual 
potentials A M = ($, A) has been commented on recently by Witten. 20 On the other hand the 
interaction of a charged particle with the field is naturally described via the potentials A M , by 
replacing the free momentum with the gauge invariant (p^ — eA^) , whereas the role of the 
{a, (3} potentials in gauge theory is not clear. 



Appendix A - Solutions of (13) in parabolic coordinates 

We solve the equations for the potentials (13). The spherical solutions yield the multipole fields, 
here we discuss a few of the parabolic solutions. Defining A, (i, 4> such that A, \i have the dimension 
of length: 21 ' 22 

xi + ix 2 — 2\f\^,e 1 ^ X3 = A — [i r — A + fx , 



then 



21,22 



v2 -ih {dxXdx+d ^ )+ ^ 

v 2 - dl + \ ee -J—(d x \a x + + 2d t ) + - d\ 



and with a = f(\)g(^)e im<t ' e ikt then 



[(d x Xd x + - 2ik) - — - — + k 2 \ + fc 2 M ] /(A) 5 (m) = . 
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Inserting a separation constant 2ikc , then 



2 

[(d x \d x - ik(l + 2c) - ^ + k 2 X] /(A) = (Al) 

2 

[d^ - ik(l -2c)-^- + feV] <?(m) = • (A2) 

The solution of (Al) is /(A) = A m/2 e- ifeA 1^1(1 + f+ c,m + l,2ik\) , where ^ is the confluent 
hypergeometric function. So the solution for a is 

a = (A^) m/2 e™* e - lfc(A+At) 1^(1 + f + c, m + 1, 2ifcA)i-Fi(l + f - c, m + 1, 2ifc/i) . (A3) 

The confluent hypergeometric function reduces to a simpler expression for integer or half-integer 
values of c . 

The case m — , c = . 

As i-Fi(l, 1,«C) = ^ we have 

_ g — ikt ^ — ik(X+fi) g2ifeAg2ife/j _ ^ — ikt^ik(X+fi) — ^ikr — ikt 

which potential substituted into (6a) yields a zero field. 
The case m — 1 , c — 1/2 . 

Inserting these values of m, c into (Al), and using the identities i-F\(2 , 2 , 2i() = e 2lC - , 
iFi(l, 2, 2tC) = ^e^ , then 

a = x /\ile i<t 'e- ik( - x+ ^ ^(2, 2, 2ifcA)iF l (l, 2, 2zfc^) e- ifet 

a = v^e^e"^^ e 2ifeA ^Me !t " e- ifet 

a =^ e ^ fcA ^Me-^ 
which is the «x of (19). 
The case m — 1, c = . 

Substituting these values for m, c, into (3.2) we have 

az = VAT* e^e-^+^iFid, 2, 2zfcA)iFi(|, 2, 2ifc^) e - lfct . 

As iFi(|,2,2i() = [Jo(C) + i«/i(C)] e 1 ^ > where J , Ji are the Bessel functions of order zero and 
one, we have 

a z = ^e^[{J +iJ l ){k\)][{J + iJ l )(k^)]e- ikt (A4) 

where we have written (Jo + iJi)(k\) as shorthand for Jo(fcA) + iJ\(k\) . For large A, vA(Jo + 
iJ\){k\) ~ \f2/irk e lkX ^ 7r / i , so that az is everywhere bounded. The potential @z such that the 
constraint (9) is satisfied is 

Pz = V^e**[(Jo - »Ji)(fcA)][(J " ^i)(M] e- lkt . (A5) 
We can determine the corresponding Debye potential ip z , such that i(d r rtpz — dtripz) = a z ■ 

^z = \ (t^-) e^JoikXjJ^kfi) + J!(fcA)J (M] • 
ik \X + n J 
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The field F z may be calculated by substituting (A4), (A5) into (8), or (A6) into (3). After 
some labour one obtains the components of Fz 



Ft 
F 2 
F 3 



-, — ikt 



(a + m) 



A[(J + iJi)(fcA) J (M] + M^o(fcA) (Jo - iJi)(kn)] \ 
iX[(J + Ui)(k\) J (M] + in[J (kX) (J - 

-2zv / AMe^[J (fcA)Ji(fcM) + Ji(kX)J (kfi)} J 



+ e 



2i</> 



-i/x[Ji(fcA)(J + iJi)(M] 
-/x[Ji(fcA)(J +iJi)(fc/x)] 



iA[(J -zJi)(A;A) Ji(kfi)} 
\[(Jo ~ iJi)(k\)] Ji(M 







(A7) 



It is interesting to see the behaviour of this field F z along the x 3 axis. For the positive x 3 axis 
p, = and A = a; 3 , and the field is just 



2(F Z ) 



ifet 



/x=0 



(J + iJ\)(kx 3 ) 
i(J + iJi){kx 3 ) 



fciC3»l 



3 -i7r/4 



7TF3 



^ikx 3 — ikt 



(A8) 



For the negative x 3 axis A = and /i = — rr 3 , and we obtain the same expression (A8) for the 
field along the negative x 3 axis. Along the x 3 axis the field Fz appears as a right-circularly 
polarized plane wave propagating in the +x 3 direction, with amplitude decaying by |a; 3 | _1 / 2 from 
the origin. Away from the x 3 axis the field decays more rapidly: the x\,x 2 plane through the 
origin is parametrized by A — p, — \j x\ + x\j2 , then substituting A — p = p into (3.7) one finds 
that the field decays by p^ 1 from the origin. Thus the field Fz is highly directional along the axis 
of propagation. 



When m — 1 , c > 1/2 , the resulting field becomes infinite at spatial infinity. 



Appendix B. Orthonormality of the basis functions £/(k) 

First we will consider the orthogonality of two left-circularly polarized waves, ?7z,(k) , £7z,(k') . Due 
to the fact that the C/i(k) , £/z,(k') are eigenstates of the angular momentum in the direction of 
propagation k or k' , and the L operator is Hermitian, then 



(£/ L (k)|£/ L (k'))=0 



when 



k^k'. 



So we need only consider the case when the waves are propagating in the same direction but with 
different frequencies, i.e. when k = {0, 0, k} and k' — {0, 0, k'} . Then recalling the Ul of (19) we 
have the inner product 



(U L (k)\U L (k')) 



= -- 4 ld*x 



,Xl+iX2, ikx3 _ £ik _ , X!+ix 2 , ikxa _ e _ ikr \ f^ £l _ iS . 2 ) e -ikx a \ 

1 - x 3 1 + x 3 J 



( X 1 -lX 2 -,r 



^3 
-ik' X3 



-ik ri 



, X\ - IX 2 \ 
1 1-X 3 ' 



-ik'x 3 „ik'r] 



e lkr ] {k( Xl + ix 2 ) e ikxs ) 



= i y d 3 x k' ((1 + x 3 ) [ e ^ k - k '^ - e i(kr-k'x 3 )^ + (i _ x 3 ) [ e ^ k - k '^ - e -«*r-«*'*a] ^ 
+ k ((1 + x 3 ) [ e ^ k - k '^ - e -i(feV-fex 3 )] + (i _ x 3 ) [ e i( - k - k '^ - e i{k ' ' r + ikx ^] ) 

= \Jd 3 x \2{k + k')e l{ - k - k '^ 

- (1 + x 3 ) [k'e l(kr - k ' X3) + ke- i{k ' r - kX3) ] - (1 - x 3 ) [k' e - ikr - ik ' X3 + ke i{k ' r+ikX3) ] 
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The first term yields the delta function (2ir) 3 k <5(k — k') . We go over to spherical coordinates 
to evaluate the remaining terms, putting £3 = cos 6 — v , the rest of the integral is 

- ^J™r 2 dr J' dv ((l + v)[k'e ir ( k - k '^+ke- ir ( k '- k ^} + (l - v) [k' e - ir ( k + k '"') 4. k e ir ( k '+ k ")j ^ 

= - — / r 2 dr / dv ((l + 1/) [fc' cos(r(k - k'v)) + k cos(r(k' - kv))] ) 
4 Jo J-i 

47T Z" 00 f 1 \ 

= — y r 2 (ir — dy) J dv i-[k' sin(r(fc — k'v)) — k sin(r(fc' — kv))] J 



i7F f r 2 dr (d k - dy ) 
Jo 



x r 2 



[cos(r(fc — k'v)) — cos(r(fc' — fci/))] 



0, 



4 

so that finally 

([/ L (k)|C/ i (k')> = (2^) 3 fc ( 5(k-k') 
with k = |k| . The other identities of (31) follow similarly. 
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